Statistical methods for assessing the effects of de novo variants on birth defects.

With the development of next-generation sequencing technology, de novo variants (DNVs) with deleterious effects can be identified and investigated for their effects on birth defects such as congenital heart disease (CHD). However, statistical power is still limited for such studies because of the small sample size due to the high cost of recruiting and sequencing samples and the low occurrence of DNVs. DNV analysis is further complicated by genetic heterogeneity across diseased individuals. Therefore, it is critical to jointly analyze DNVs with other types of genomic/biological information to improve statistical power to identify genes associated with birth defects. In this review, we discuss the general workflow, recent developments in statistical methods, and future directions for DNV analysis.

Clinical Biochemistry of Serum Troponin.

Accurate measurement and interpretation of serum levels of troponin (Tn) is a central part of the clinical workup of a patient presenting with chest pain suspicious for acute coronary syndrome (ACS). Knowledge of the molecular characteristics of the troponin complex and test characteristics of troponin measurement assays allows for a deeper understanding of causes of false positive and false negative test results in myocardial injury. In this review, we discuss the molecular structure and functions of the constituent proteins of the troponin complex (TnT, TnC, and TnI); review the different isoforms of Tn and where they are from; survey the evolution of clinical Tn assays, ranging from first-generation to high-sensitivity (hs); provide a primer on statistical interpretation of assay results based on different clinical settings; and discuss potential causes of false results. We also summarize the advances in technologies that may lead to the development of future Tn assays, including the development of point of care assays and wearable Tn sensors for real-time continuous measurement.

Mutation of key signaling regulators of cerebrovascular development in vein of Galen malformations.

To elucidate the pathogenesis of vein of Galen malformations (VOGMs), the most common and most severe of congenital brain arteriovenous malformations, we performed an integrated analysis of 310 VOGM proband-family exomes and 336,326 human cerebrovasculature single-cell transcriptomes. We found the Ras suppressor p120 RasGAP (RASA1) harbored a genome-wide significant burden of loss-of-function de novo variants (2042.5-fold, p = 4.79 x 10-7). Rare, damaging transmitted variants were enriched in Ephrin receptor-B4 (EPHB4) (17.5-fold, p = 1.22 x 10-5), which cooperates with p120 RasGAP to regulate vascular development. Additional probands had damaging variants in ACVRL1, NOTCH1, ITGB1, and PTPN11. ACVRL1 variants were also identified in a multi-generational VOGM pedigree. Integrative genomic analysis defined developing endothelial cells as a likely spatio-temporal locus of VOGM pathophysiology. Mice expressing a VOGM-specific EPHB4 kinase-domain missense variant (Phe867Leu) exhibited disrupted developmental angiogenesis and impaired hierarchical development of arterial-capillary-venous networks, but only in the presence of a "second-hit" allele. These results illuminate human arterio-venous development and VOGM pathobiology and have implications for patients and their families.

Nonexercise machine learning models for maximal oxygen uptake prediction in national population surveys.

OBJECTIVE: Nonexercise algorithms are cost-effective methods to estimate cardiorespiratory fitness (CRF), but the existing models have limitations in generalizability and predictive power. This study aims to improve the nonexercise algorithms using machine learning (ML) methods and data from US national population surveys. MATERIALS AND METHODS: We used the 1999-2004 data from the National Health and Nutrition Examination Survey (NHANES). Maximal oxygen uptake (VO2 max), measured through a submaximal exercise test, served as the gold standard measure for CRF in this study. We applied multiple ML algorithms to build 2 models: a parsimonious model using commonly available interview and examination data, and an extended model additionally incorporating variables from Dual-Energy X-ray Absorptiometry (DEXA) and standard laboratory tests in clinical practice. Key predictors were identified using Shapley additive explanation (SHAP). RESULTS: Among the 5668 NHANES participants in the study population, 49.9% were women and the mean (SD) age was 32.5 years (10.0). The light gradient boosting machine (LightGBM) had the best performance across multiple types of supervised ML algorithms. Compared with the best existing nonexercise algorithms that could be applied to the NHANES, the parsimonious LightGBM model (RMSE: 8.51 ml/kg/min [95% CI: 7.73-9.33]) and the extended LightGBM model (RMSE: 8.26 ml/kg/min [95% CI: 7.44-9.09]) significantly reduced the error by 15% and 12% (P < .001 for both), respectively. DISCUSSION: The integration of ML and national data source presents a novel approach for estimating cardiovascular fitness. This method provides valuable insights for cardiovascular disease risk classification and clinical decision-making, ultimately leading to improved health outcomes. CONCLUSION: Our nonexercise models provide improved accuracy in estimating VO2 max within NHANES data as compared to existing nonexercise algorithms.

Genetic dysregulation of an endothelial Ras signaling network in vein of Galen malformations.

To elucidate the pathogenesis of vein of Galen malformations (VOGMs), the most common and severe congenital brain arteriovenous malformation, we performed an integrated analysis of 310 VOGM proband-family exomes and 336,326 human cerebrovasculature single-cell transcriptomes. We found the Ras suppressor p120 RasGAP ( RASA1 ) harbored a genome-wide significant burden of loss-of-function de novo variants (p=4.79×10 -7 ). Rare, damaging transmitted variants were enriched in Ephrin receptor-B4 ( EPHB4 ) (p=1.22×10 -5 ), which cooperates with p120 RasGAP to limit Ras activation. Other probands had pathogenic variants in ACVRL1 , NOTCH1 , ITGB1 , and PTPN11 . ACVRL1 variants were also identified in a multi-generational VOGM pedigree. Integrative genomics defined developing endothelial cells as a key spatio-temporal locus of VOGM pathophysiology. Mice expressing a VOGM-specific EPHB4 kinase-domain missense variant exhibited constitutive endothelial Ras/ERK/MAPK activation and impaired hierarchical development of angiogenesis-regulated arterial-capillary-venous networks, but only when carrying a "second-hit" allele. These results illuminate human arterio-venous development and VOGM pathobiology and have clinical implications.

Uterine leiomyosarcomas harboring MAP2K4 gene amplification are sensitive in vivo to PLX8725, a novel MAP2K4 inhibitor.

INTRODUCTION: Uterine leiomyosarcomas (uLMS) are rare, highly aggressive tumors. Up to 30% of uLMS may harbor gain of function (GOF) in the MAP2K4 gene, important for tumor cell proliferation, differentiation and metastasis. We investigated the in vivo activity of a novel MAP2K4 inhibitor, PLX8725, against uLMS harboring MAP2K4 gene-amplification. METHODS: Two fully characterized uLMS (i.e., LEY-11 and LEY-16) were grafted into female CB-17/SCID mice. Treatments with control vehicle or PLX8725 (50 mg/kg) were given via oral gavage daily on weekdays for up to 60 days. Tumor volume differences were calculated with two-way ANOVA. Pharmacokinetic (PK) and mechanistic studies of PLX8725 in uLMS PDX models were also performed. RESULTS: Both uLMS tumors evaluated demonstrated GOF in MAP2K4 (i.e., 3 CNV in both LEY-11 and LEY-16). Tumor growth inhibition was significantly greater in both PDX LEY-11 and PDX LEY-16 treated with PLX8725 when compared to controls (p < 0.001). Median overall survival was also significantly longer in both PDX LEY-11 (p = 0.0047) and PDX LEY-16 (p = 0.0058) treatment cohorts when compared to controls. PLX8725 oral treatment was well tolerated, and PK studies demonstrated that oral PLX8725 gives extended exposure in mice. Ex vivo tumor samples after PLX8725 exposure decreased phosphorylated-ATR, JNK and p38, and increased expression of apoptotic molecules on western blot. CONCLUSION: PLX8725 demonstrates promising in vivo activity against PDX models of uLMS harboring GOF alterations in the MAP2K4 gene with tolerable toxicity. Phase I trials of PLX8725 in advanced, recurrent, chemotherapy-resistant uLMS patients are warranted.

Multiomic analyses implicate a neurodevelopmental program in the pathogenesis of cerebral arachnoid cysts.

Cerebral arachnoid cysts (ACs) are one of the most common and poorly understood types of developmental brain lesion. To begin to elucidate AC pathogenesis, we performed an integrated analysis of 617 patient-parent (trio) exomes, 152,898 human brain and mouse meningeal single-cell RNA sequencing transcriptomes and natural language processing data of patient medical records. We found that damaging de novo variants (DNVs) were highly enriched in patients with ACs compared with healthy individuals (P = 1.57 × 10-33). Seven genes harbored an exome-wide significant DNV burden. AC-associated genes were enriched for chromatin modifiers and converged in midgestational transcription networks essential for neural and meningeal development. Unsupervised clustering of patient phenotypes identified four AC subtypes and clinical severity correlated with the presence of a damaging DNV. These data provide insights into the coordinated regulation of brain and meningeal development and implicate epigenomic dysregulation due to DNVs in AC pathogenesis. Our results provide a preliminary indication that, in the appropriate clinical context, ACs may be considered radiographic harbingers of neurodevelopmental pathology warranting genetic testing and neurobehavioral follow-up. These data highlight the utility of a systems-level, multiomics approach to elucidate sporadic structural brain disease.

Elimusertib (BAY1895344), a novel ATR inhibitor, demonstrates in vivo activity in ATRX mutated models of uterine leiomyosarcoma.

INTRODUCTION: Uterine leiomyosarcoma (uLMS) is a rare, highly aggressive malignancy. Recent data suggest 50% of uLMS may harbor alterations in the ATRX gene and such mutations may confer sensitivity to ataxia-telangiectasia-and-Rad3-related (ATR) kinase inhibitors. We sought to investigate the in vivo activity of Elimusertib (BAY1895344), a novel ATR-inhibitor, against ATRX-mutated uLMS patient-derived xenografts (PDXs). METHODS: Two fully characterized uLMS (i.e., LEY-11 and LEY-16) were grafted into female CB-17/SCID mice. Treatments with control vehicle or BAY1895344 (20 mg/kg dosed twice daily 3 days on 4 days off) were given via oral gavage and tumor measurements as well as weights obtained twice weekly. Tumor volume differences were calculated with a two-way ANOVA. Mechanistic studies were performed ex vivo using BAY1895344 treated uLMS tumor samples by western blot analysis. RESULTS: Both PDX LEY-11 and PDX LEY-16 harboring ATRX gene mutations demonstrated an aggressive behavior in vivo (i.e., control mice were euthanized on average at day 12.5 for PDX LEY-11 and at day 33 for PDX LEY-16). In both tumor models BAY1895344 20 mg/kg dosed with an intermittent oral schedule was able to induce significant growth inhibition compared to vehicle control treatment (p < 0.001 for both LEY-11 and LEY-16) and prolong median overall survival [PDX LEY-11 (12.5 vs. 42 days, p < 0.001) and PDX LEY-16 (33 vs. 60 days, p < 0.001)]. There were not significant changes in weight between treatment and controls. By western blot assays BAY1895344 exposure decreased phosphorylated-ATR and increased expression of apoptotic molecules in LMS PDXs. CONCLUSIONS: BAY1895344 demonstrates promising in vivo activity against biologically aggressive PDX models of uLMS harboring ATRX mutations, with no significant toxicity. Clinical trials of BAY1895344 in uLMS patients are warranted.

A Novel Missense Mutation in ERCC8 Co-Segregates with Cerebellar Ataxia in a Consanguineous Pakistani Family.

Autosomal-recessive cerebellar ataxias (ARCAs) are heterogeneous rare disorders mainly affecting the cerebellum and manifest as movement disorders in children and young adults. To date, ARCA causing mutations have been identified in nearly 100 genes; however, they account for less than 50% of all cases. We studied a multiplex, consanguineous Pakistani family presenting with a slowly progressive gait ataxia, body imbalance, and dysarthria. Cerebellar atrophy was identified by magnetic resonance imaging of brain. Using whole exome sequencing, a novel homozygous missense mutation ERCC8:c.176T>C (p.M59T) was identified that co-segregated with the disease. Previous studies have identified homozygous mutations in ERCC8 as causal for Cockayne Syndrome type A (CSA), a UV light-sensitive syndrome, and several ARCAs. ERCC8 plays critical roles in the nucleotide excision repair complex. The p.M59T, a substitution mutation, is located in a highly conserved WD1 beta-transducin repeat motif. In silico modeling showed that the structure of this protein is significantly affected by the p.M59T mutation, likely impairing complex formation and protein-protein interactions. In cultured cells, the p.M59T mutation significantly lowered protein stability compared to wildtype ERCC8 protein. These findings expand the role of ERCC8 mutations in ARCAs and indicate that ERCC8-related mutations should be considered in the differential diagnosis of ARCAs.

Network assisted analysis of de novo variants using protein-protein interaction information identified 46 candidate genes for congenital heart disease.

De novo variants (DNVs) with deleterious effects have proved informative in identifying risk genes for early-onset diseases such as congenital heart disease (CHD). A number of statistical methods have been proposed for family-based studies or case/control studies to identify risk genes by screening genes with more DNVs than expected by chance in Whole Exome Sequencing (WES) studies. However, the statistical power is still limited for cohorts with thousands of subjects. Under the hypothesis that connected genes in protein-protein interaction (PPI) networks are more likely to share similar disease association status, we developed a Markov Random Field model that can leverage information from publicly available PPI databases to increase power in identifying risk genes. We identified 46 candidate genes with at least 1 DNV in the CHD study cohort, including 18 known human CHD genes and 35 highly expressed genes in mouse developing heart. Our results may shed new insight on the shared protein functionality among risk genes for CHD.

Whole-exome sequencing reveals damaging gene variants associated with hypoalphalipoproteinemia.

Low levels of high density lipoprotein-cholesterol (HDL-C) are associated with an elevated risk of arteriosclerotic coronary heart disease. Heritability of HDL-C levels is high. In this research discovery study, we used whole-exome sequencing to identify damaging gene variants that may play significant roles in determining HDL-C levels. We studied 204 individuals with a mean HDL-C level of 27.8 ± 6.4 mg/dl (range: 4-36 mg/dl). Data were analyzed by statistical gene burden testing and by filtering against candidate gene lists. We found 120 occurrences of probably damaging variants (116 heterozygous; four homozygous) among 45 of 104 recognized HDL candidate genes. Those with the highest prevalence of damaging variants were ABCA1 (n = 20), STAB1 (n = 9), OSBPL1A (n = 8), CPS1 (n = 8), CD36 (n = 7), LRP1 (n = 6), ABCA8 (n = 6), GOT2 (n = 5), AMPD3 (n = 5), WWOX (n = 4), and IRS1 (n = 4). Binomial analysis for damaging missense or loss-of-function variants identified the ABCA1 and LDLR genes at genome-wide significance. In conclusion, whole-exome sequencing of individuals with low HDL-C showed the burden of damaging rare variants in the ABCA1 and LDLR genes is particularly high and revealed numerous occurrences in HDL candidate genes, including many genes identified in genome-wide association study reports. Many of these genes are involved in cancer biology, which accords with epidemiologic findings of the association of HDL deficiency with increased risk of cancer, thus presenting a new area of interest in HDL genomics.

Impaired neurogenesis alters brain biomechanics in a neuroprogenitor-based genetic subtype of congenital hydrocephalus.

Hydrocephalus, characterized by cerebral ventricular dilatation, is routinely attributed to primary defects in cerebrospinal fluid (CSF) homeostasis. This fosters CSF shunting as the leading reason for brain surgery in children despite considerable disease heterogeneity. In this study, by integrating human brain transcriptomics with whole-exome sequencing of 483 patients with congenital hydrocephalus (CH), we found convergence of CH risk genes in embryonic neuroepithelial stem cells. Of all CH risk genes, TRIM71/lin-41 harbors the most de novo mutations and is most specifically expressed in neuroepithelial cells. Mice harboring neuroepithelial cell-specific Trim71 deletion or CH-specific Trim71 mutation exhibit prenatal hydrocephalus. CH mutations disrupt TRIM71 binding to its RNA targets, causing premature neuroepithelial cell differentiation and reduced neurogenesis. Cortical hypoplasia leads to a hypercompliant cortex and secondary ventricular enlargement without primary defects in CSF circulation. These data highlight the importance of precisely regulated neuroepithelial cell fate for normal brain-CSF biomechanics and support a clinically relevant neuroprogenitor-based paradigm of CH.

Mutation spectrum of congenital heart disease in a consanguineous Turkish population.

BACKGROUNDS: While many studies agree that consanguinity increases the rate of congenital heart disease (CHD), few genome analyses have been conducted with consanguineous CHD cohorts. METHODS: We recruited 73 CHD probands from consanguineous families in Turkey and used whole-exome sequencing (WES) to identify genetic lesions in these patients. RESULTS: On average, each patient had 6.95 rare damaging homozygous variants, 0.68 of which are loss-of-function (LoF) variants. Seven patients (9.6%) carried damaging homozygous variants in five causal CHD genes. Six of those patients exhibited laterality defects (six HTX and one D-TGA). Three additional patients (4.1%) harbored other types of CHD-associated genomic alterations, which overall explained 13.7% (10/73) of the cohort. The contribution from recessive variants in our cohort is higher than 1.8% reported from a cohort of 2871 CHD subjects where 5.6% of subjects met the criteria for consanguinity. CONCLUSIONS: Our WES screen of a Turkish consanguineous population with structural CHD revealed its unique genetic architecture. Six of seven damaging homozygous variants in CHD causal genes occur in the setting of laterality defects implies a strong contribution from consanguinity to these defects specifically. Our study thus provided valuable information about the genetic landscape of CHD in consanguineous families in Turkey.

A phase 2 evaluation of pembrolizumab for recurrent Lynch-like versus sporadic endometrial cancers with microsatellite instability.

BACKGROUND: Microsatellite instability-high (MSI-H)/mismatch repair deficiency (dMMR) is a biomarker for responses to immune checkpoint inhibitors (ICIs). Whether mechanisms underlying microsatellite instability alter responses to ICIs is unclear. This article reports data from a prospective phase 2 pilot study of pembrolizumab in patients with recurrent MSI-H endometrial cancer (EC) analyzed by whole exome sequencing (WES) and potential mechanisms of primary/secondary ICI resistance (NCT02899793). METHODS: Patients with measurable MSI-H/dMMR EC confirmed by polymerase chain reaction/immunohistochemistry were evaluated by WES and received 200 mg of pembrolizumab every 3 weeks for ≤2 years. The primary end point was the objective response rate (ORR). Secondary end points included progression-free survival (PFS) and overall survival (OS). RESULTS: Twenty-five patients (24 evaluable) were treated. Six patients (25%) harbored Lynch/Lynch-like tumors, whereas 18 (75%) had sporadic EC. The tumor mutation burden was higher in Lynch-like tumors (median, 2939 mutations/megabase [Mut/Mb]; interquartile range [IQR], 867-5108 Mut/Mb) than sporadic tumors (median, 604 Mut/Mb; IQR, 411-798 Mut/Mb; P = .0076). The ORR was 100% in Lynch/Lynch-like patients but only 44% in sporadic patients (P = .024). The 3-year PFS and OS proportions were 100% versus 30% (P = .017) and 100% versus 43% (P = .043), respectively. CONCLUSIONS: This study suggests prognostic significance of Lynch-like cancers versus sporadic MSI-H/dMMR ECs for ORR, PFS, and OS when patients are treated with pembrolizumab. Larger confirmatory studies in ECs and other MSI-H/dMMR tumors are necessary. Defective antigen processing/presentation and deranged induction in interferon responses serve as mechanisms of resistance in sporadic MSI-H ECs. Oligoprogression in MSI-H/dMMR patients appears salvageable with surgical resection and/or local treatment and the continuation of pembrolizumab off study. Clinical studies evaluating separate MSI-H/dMMR EC subtypes treated with ICIs are warranted.

PPIL4 is essential for brain angiogenesis and implicated in intracranial aneurysms in humans.

Intracranial aneurysm (IA) rupture leads to subarachnoid hemorrhage, a sudden-onset disease that often causes death or severe disability. Although genome-wide association studies have identified common genetic variants that increase IA risk moderately, the contribution of variants with large effect remains poorly defined. Using whole-exome sequencing, we identified significant enrichment of rare, deleterious mutations in PPIL4, encoding peptidyl-prolyl cis-trans isomerase-like 4, in both familial and index IA cases. Ppil4 depletion in vertebrate models causes intracerebral hemorrhage, defects in cerebrovascular morphology and impaired Wnt signaling. Wild-type, but not IA-mutant, PPIL4 potentiates Wnt signaling by binding JMJD6, a known angiogenesis regulator and Wnt activator. These findings identify a novel PPIL4-dependent Wnt signaling mechanism involved in brain-specific angiogenesis and maintenance of cerebrovascular integrity and implicate PPIL4 gene mutations in the pathogenesis of IA.

M-DATA: A statistical approach to jointly analyzing de novo mutations for multiple traits.

Recent studies have demonstrated that multiple early-onset diseases have shared risk genes, based on findings from de novo mutations (DNMs). Therefore, we may leverage information from one trait to improve statistical power to identify genes for another trait. However, there are few methods that can jointly analyze DNMs from multiple traits. In this study, we develop a framework called M-DATA (Multi-trait framework for De novo mutation Association Test with Annotations) to increase the statistical power of association analysis by integrating data from multiple correlated traits and their functional annotations. Using the number of DNMs from multiple diseases, we develop a method based on an Expectation-Maximization algorithm to both infer the degree of association between two diseases as well as to estimate the gene association probability for each disease. We apply our method to a case study of jointly analyzing data from congenital heart disease (CHD) and autism. Our method was able to identify 23 genes for CHD from joint analysis, including 12 novel genes, which is substantially more than single-trait analysis, leading to novel insights into CHD disease etiology.

Molecular Genetics and Complex Inheritance of Congenital Heart Disease.

Congenital heart disease (CHD) is the most common congenital malformation and the leading cause of mortality therein. Genetic etiologies contribute to an estimated 90% of CHD cases, but so far, a molecular diagnosis remains unsolved in up to 55% of patients. Copy number variations and aneuploidy account for ~23% of cases overall, and high-throughput genomic technologies have revealed additional types of genetic variation in CHD. The first CHD risk genotypes identified through high-throughput sequencing were de novo mutations, many of which occur in chromatin modifying genes. Murine models of cardiogenesis further support the damaging nature of chromatin modifying CHD mutations. Transmitted mutations have also been identified through sequencing of population scale CHD cohorts, and many transmitted mutations are enriched in cilia genes and Notch or VEGF pathway genes. While we have come a long way in identifying the causes of CHD, more work is required to end the diagnostic odyssey for all CHD families. Complex genetic explanations of CHD are emerging but will require increasingly sophisticated analysis strategies applied to very large CHD cohorts before they can come to fruition in providing molecular diagnoses to genetically unsolved patients. In this review, we discuss the genetic architecture of CHD and biological pathways involved in its pathogenesis.

DIAPH1 Variants in Non-East Asian Patients With Sporadic Moyamoya Disease.

Importance: Moyamoya disease (MMD), a progressive vasculopathy leading to narrowing and ultimate occlusion of the intracranial internal carotid arteries, is a cause of childhood stroke. The cause of MMD is poorly understood, but genetic factors play a role. Several familial forms of MMD have been identified, but the cause of most cases remains elusive, especially among non-East Asian individuals. Objective: To assess whether ultrarare de novo and rare, damaging transmitted variants with large effect sizes are associated with MMD risk. Design, Setting, and Participants: A genetic association study was conducted using whole-exome sequencing case-parent MMD trios in a small discovery cohort collected over 3.5 years (2016-2019); data were analyzed in 2020. Medical records from US hospitals spanning a range of 1 month to 1.5 years were reviewed for phenotyping. Exomes from a larger validation cohort were analyzed to identify additional rare, large-effect variants in the top candidate gene. Participants included patients with MMD and, when available, their parents. All participants who met criteria and were presented with the option to join the study agreed to do so; none were excluded. Twenty-four probands (22 trios and 2 singletons) composed the discovery cohort, and 84 probands (29 trios and 55 singletons) composed the validation cohort. Main Outcomes and Measures: Gene variants were identified and filtered using stringent criteria. Enrichment and case-control tests assessed gene-level variant burden. In silico modeling estimated the probability of variant association with protein structure. Integrative genomics assessed expression patterns of MMD risk genes derived from single-cell RNA sequencing data of human and mouse brain tissue. Results: Of the 24 patients in the discovery cohort, 14 (58.3%) were men and 18 (75.0%) were of European ancestry. Three of 24 discovery cohort probands contained 2 do novo (1-tailed Poisson P = 1.1 × 10-6) and 1 rare, transmitted damaging variant (12.5% of cases) in DIAPH1 (mammalian diaphanous-1), a key regulator of actin remodeling in vascular cells and platelets. Four additional ultrarare damaging heterozygous DIAPH1 variants (3 unphased) were identified in 3 other patients in an 84-proband validation cohort (73.8% female, 77.4% European). All 6 patients were non-East Asian. Compound heterozygous variants were identified in ena/vasodilator-stimulated phosphoproteinlike protein EVL, a mammalian diaphanous-1 interactor that regulates actin polymerization. DIAPH1 and EVL mutant probands had severe, bilateral MMD associated with transfusion-dependent thrombocytopenia. DIAPH1 and other MMD risk genes are enriched in mural cells of midgestational human brain. The DIAPH1 coexpression network converges in vascular cell actin cytoskeleton regulatory pathways. Conclusions and Relevance: These findings provide the largest collection to date of non-East Asian individuals with sporadic MMD harboring pathogenic variants in the same gene. The results suggest that DIAPH1 is a novel MMD risk gene and impaired vascular cell actin remodeling in MMD pathogenesis, with diagnostic and therapeutic ramifications.